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Использование особенностей возмущений 
сингулярных чисел матрицы цифрового 
изображения для обнаружения 

его фальсификации 


Предлагается новый подход к решению проблемы обнаружения фальсификации цифрового 
изображения, основанный на анализе возмущений сингулярных чисел блоков матрицы изображения. 
Новый подход дает возможность локализовать область несанкционированного изменения исходного 
изображения. 


Введение 


Доказательство подлинности цифровых изображений (ЦИ), обнаружение их 
фальсификации является важной задачей во многих областях человеческой 
деятельности (работа с рентгенографическими и томографическими изображениями 
в медицине, судебные разбирательства, пресса, научные публикации и т.д.). 
Возрастающая актуальность рассматриваемой задачи в настоящий момент в 
значительной степени обязана последним успехам в технике синтеза изображений, а 
также общедоступности редактирующего изображения программного обеспечения. 

Одним из эффективных инструментов доказательства подлинности ЦИ 
является использование техники цифровых водяных знаков (ЦВЗ) [1-3]. Основной 
недостаток этого подхода состоит в том, что ЦВЗ должен быть погружен в цифровой 
объект во время создания этого объекта, что ограничивает область применения 
метода только для механизмов генерации ЦИ, имеющих встроенные возможности 
погружения ЦВЗ. 

В последнее время были предложены новые подходы к решению рас- 
сматриваемой задачи: метод автоматического детектирования характерных для 
фальсификации статистических корреляций, являющихся результатом передискре- 
тизации (те-затрИп?) изображений [4]; техника, основанная на оценке место- 
положения источника света, учитывающая световые различия, которые имеют место 
в случае, если изображение является комбинацией нескольких [5|; метод, основанный на 
выявлении результатов элементарных операций обработки изображения, используемых 
при его фальсификации (масштабирование, поворот, сглаживание и т.д.) [6]. 


Постановка задачи и цель исследования 


Каждый из указанных выше методов имеет конкретную область применения. 
Большинство из методов не гарантируют обнаружение фальсификации ЦИ при ее 
наличии, а при обнаружении — не локализуют ее область. Это заставляет искать 
новые математические инструменты и подходы к решению рассматриваемой задачи 
в целом. 


«Штучний 1нтелект» 12008 145 


| Кобозева А.А. 


Глобальной целью автора является создание единого математического 
подхода к решению проблемы доказательства подлинности ЦИ, детектирования его 
фальсификации, независимо от вида фальсифицирующего воздействия. 

Любая фальсификация ЦИ представляет из себя возмущение [7] матрицы 
(матриц) изображения, а значит, может быть представлена как совокупность 
возмущений сингулярных чисел (СЧ) и сингулярных векторов (СВ) исходной 
матрицы, однозначно определяемых ее нормальным сингулярным разложением [8]. 
Глобальная задача заключается в том, чтобы установить характерные признаки 
этих возмущений (или свойства СЧ и СВ), наличие или отсутствие которых даст 
возможность не только отделить фальсифицированное изображение от ориги- 
нального, но и локализовать область фальсификации. 

Целью настоящей работы является создание начальных основ нового 
математического подхода к решению рассматриваемой задачи, базирующегося на 
теории возмущений и матричном анализе. 


Характерные особенности сингулярных чисел 


Пусть исходное ЦИ, хранящееся в некотором формате без потерь, например в 
формате ТТЕ, разбивается стандартным образом на блоки 8х8 [9]. Если для каждого 
блока изображения определить сингулярный спектр (множество всех СЧ), то 
оказывается, что менее 10 % от общего числа блоков имеют нулевые СЧ (табл. 1). 


Таблица 1 — Результаты ЗУ блоков 8х8 изображений, хранимых без потерь 


Общее | Количество | Количество | Количество | Количество Кол-во 
число блоков, блоков, блоков, блоков, блоков, 
Изображение блоков | имеющих 1 | имеющих 2 | имеющих 3 | имеющих 4 имеющих 
(ОЧБ) | нулевое СЧ нулевых нулевых нулевых | нулевые СЧ, 
СЧ СЧ СЧ по 
отношению к 
ОЧБ (в %) 
РОЧТ.ТЕ 1080 10 0 0 0 0,9 
САМЕКАМАМ.ТТЕ 1024 52 33 1 0 8,4 
ТТВЕ. ТЕ 725 9 0 1 0 1,4 
МООМ.ТТЕ 2948 153 51 34 14 9 
СЕБТ. ТЕ 437 0 0 0 0 0 


Данный факт не случаен. Ранг любой матрицы определяется количеством ее 
ненулевых СЧ [7], а значит, наличие нулей в сингулярном спектре будет говорить о 
том, что число ее линейно независимых строк (столбцов) меньше размерности. 
Однако для произвольного изображения, даже с учетом коррелированности значений 
яркости пикселей, вероятность того, что строки (столбцы) очередного блока окажутся 
линейно зависимыми, невелика. Чаще всего это возникает в случае коллинеарности 
двух векторов, либо когда рассматриваемый блок соответствует очень темной области 
изображения без контуров (все элементы матрицы блока близки к 0). 

В настоящее время большинство цифровых камер используют для хранения 
изображений схему ЛРЕС (с потерями), основанную на дискретном косинусном 
преобразовании (ОСТ), включающую 3 основных шага: 

1. ОСТ. Матрица ЦИ разбивается на блоки 8х8, для каждого из которых 

выполняется сдвиг диапазона значений на 128 и прямое ОСТ. 

2. Квантование. ОСТ-коэффициенты квантуются при помощи массива норма- 

лизации и округляются до ближайшего целого. 

3. Энтропийное кодирование. Наиболее часто здесь используется процедура 

кодирования Хаффмана [9]. 
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Восстановление изображения включает в себя шаги, обратные к пере- 
численным выше, в обратном порядке. 

Поскольку шаг 3 РЕС-кодирования является обратимым, он не представляет для 
нас интереса с точки зрения определения характерных особенностей СЧ блоков 
восстановленного изображения, равно как и сдвиг диапазона значений яркости на 128. 
Поэтому в дальнейшем энтропийное кодирование и декодирование отдельно не 
рассматриваются, происходящий сдвиг диапазона для простоты изложения опускается. 

Квантование коэффициентов ОСТ является необратимой процедурой и при- 
водит к некоторым особенностям СЧ блоков. 

Пусть исходное изображение подверглось ЛРЕС-сжатию. Проведем для него 
операцию частичного восстановления (ЧВ), которая включает в себя: 1) энтропийное 
декодирование; 2) умножение полученных после первого шага коэффициентов на 
соответствующие элементы массива нормализации; 3) применение к полученным 
после второго шага коэффициентам обратного ОСТ, но без последующего округления. 
У полученной матрицы практически все блоки содержат нулевые СЧ, причем таких 
значений в блоках будет достаточно много (табл. 2 - в колонках 3 -— 11 
располагаются количества блоков матрицы, имеющих т нулевых СЧ). Такая 
ситуация закономерна. После квантования и округления коэффициентов ОСТ блоков 
многие из них, отвечающие высоким и средним частотам, обнулятся, оставаясь 
нулями после ЧВ. Существует соответствие между элементами дискретного 


преобразования Фурье и сингулярными тройками (©,,и,,у, } матрицы изображения, 


где о,и‚,у, — СЧ и отвечающие ему левый и правый СВ соответственно [7]. 


Сингулярные тройки наибольших СЧ максимально отвечают низкочастотным 
составляющим сигнала-изображения. По мере уменьшения СЧ происходит 
увеличение доли средних и высоких частот, а вклад низких становится все меньше. 
Наименьшие СЧ отвечают высокочастотным составляющим двумерного цифрового 


А 
сигнала [8]. Если Ё, = Усил — аппроксимация ранга К изображения с матрицей РЕ, 
: 1=1 
В = Усил — дополнением к аппроксимации Е,, 5, =си,у, — К-ая состав- 
1-1 
ляющая Е, то дополнения к аппроксимации и составляющие изображения при 
достаточно большом А практически не содержат низкочастотных составляющих [8]. 
Поскольку после ЧВ подавляющее большинство коэффициентов ОСТ, отвечающих 
высокочастотным (и средним) составляющим, становятся нулевыми, то нулевыми 
будут наименьшие (а возможно, и средние по величине) СЧ матриц блоков. 


Таблица 2 — Результаты сингулярного разложения блоков 8х8 частично вос- 
становленных изображений 


Общее т Кол-во блоков, 
число у которых нулевых 
Изображение блоков | 8 7 6 5 4 3 211101 СЧ больше двух, 
по отношению к 
ОЧБ (в %) 
РОЧТ./ТЕ 1080 0 473 | 351 | 183 | 64 9 01010 100 
САМЕВКАМАМ.ТШЕ | 1024 0 472 | 69 | 72 | 123 | 167 | 99 | 20 |2 88 
ИВЕ ТЕ 725 20 | 168 | 112 | 117 | 141 | 123 | 431110 94 
МООМ.ТТЕ 2948 |218 | 1624 | 464 | 387 | 182 | 68151010 99,8 
СЕШС/ТТЕ 437 0 303 | 31| 17 | 54 | 28141000 99 


«Штучний 1нтелект» 12008 147 


| Кобозева А.А. 


Заметим, что, чем меньше нулевых СЧ в рассматриваемом блоке, тем больше 
линий контура он содержит. Действительно, наличие контуров в блоке говорит о 
значительной высокочастотной составляющей в сигнале, отвечающем этому блоку. 
Тогда коэффициенты ОСТ, соответствующие высоким и средним частотам, будут 
сравнительно большими и могут остаться ненулевыми после квантования и ЧВ, а 
значит, внесут свой вклад не только в максимальные СЧ. 

Для наглядного представления справедливости вышесказанного рассмотрим 
изображение САМЕВАМАМ.ТТЕ (рис. [(а)), подвергшееся 7РЕС-сжатию и ЧВ. На 
рис. 1(6) представлена матрица, каждый элемент которой определяет количество 
нулевых СЧ в соответствующем блоке (далее такую матрицу будем называть 
матрицей нулевых СЧ блоков (МНСЧБ)) ЧВ-изображения. Элементы, значения 
которых 0, 1, 2, 3, выделены более насыщенным тоном, что позволяет наглядно 
увидеть соответствие между контурами исходного изображения и блоками, содер- 
жащими малое количество нулевых СЧ. 


Рисунок 1 — Исходное изображение (а) и МНСЧБ после частичного 
восстановления (6) 


Таблица 3 — Результаты сингулярного разложения блоков 8х8 полностью 
восстановленных изображений 


т Кол-во блоков, 
Общее у которых 
Изображение число 8 7 6 5 4 3 2 1 0 нулевых СЧ 
блоков больше двух, по 
отношению к 
ОЧБ (в %) 
РОЧТ/Е 1080 0 473 | 126 | 125 | 110 | 93 | 93 | 29 | 31 86 
САМЕКАМАМ.ТТЕ | _1024 0 | 472 [51 | 37 | 58 | 70 [109 | 81 | 146 67 
ИВЕ ТЕ 725 20 | 168 | 47 | 42 | 59 | 55 | 80 [32 | 222 54 
МООМ. ТЕ 2948 _| 218 | 1624 | 250 | 191 | 195 | 165 | 138 | 63 | 104 90 
СЕР. ТЕ 437 0 303_| 19 7 7 14 | 25 [17| 45 80 


Пусть исходное изображение, подвергшееся ЛРЕС-сжатию, восстанавливается 
полностью. Это означает, что после ЧВ все значения яркости пикселей округляются 
до целых и вводятся в диапазон [0,255]. Это действие возмутит матрицу 
изображения, полученную после ЧВ, определенным образом изменится количество 
нулевых СЧ в блоках (табл. 3). Там, где после частичного восстановления не было 
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элементов, значительно меньших 0 или больших 255, возмущение матрицы будет 
небольшим. В соответствии с соотношением 


<|Е 


| -о' в 1=ЪП, 


имеющим место для произвольной матрицы размерности пхи [7], где о,, о: — СЧ 


М, - 
спектральная матричная норма, СЧ являются нечувствительными к возмущающим 
воздействиям. Нулевые СЧ блоков матрицы ЧВ-изображения хоть и станут не 
нулями после полного восстановления, но их значения будут сравнимы с погрешностью 
округления, что не характерно для блоков исходного изображения. Для иллюстрации 
рассмотрим блок (27,8) изображения МООМ.ТЕ. В исходном изображении син- 
гулярный спектр этого блока имеет вид: 

950.9198, 164.1247, 61.7382, 17.0963, 7.1338, 4.0991, 1.5753, 0.6998. 
После ГРЕС-сжатия и частичного восстановления СЧ блока равны: 

954.1724, 169.0268, 86.5101, 0.0000, 0.0000, 0.0000, 0.0000, 0.0000, 
а после полного восстановления: 

947.8175, 167.5342 , 89.1161, 1.6238, 0.7192, 0.5034, 0.1063, 0.0160. 


исходной И возмущенной матриц соответственно, Е - матрица возмущений, 


5 
Рисунок 2 — МНСЧЬБ, отвечающая исходному изображению (а) и полностью 
восстановленному изображению (6) 


Хотя последние 5 СЧ стали ненулевыми, но их значения незначительно отличаются 
от нуля и друг от друга и явно отличаются от трех максимальных. Такая особенность дает 
возможность различать блоки изображения, полностью восстановленного после Л/РЕС- 
сжатия, и блоки исходного изображения. 

Наиболее заметным различие между совокупным исходным изображением и 
полностью восстановленным после ЛРЕС-сжатия будет при сравнении их МНСЧБ. 
Типичная картина представлена на рис. 2. (изображение САМЕКАМАМ.Т!®). Хотя 
значения 0, 1, 2, 3 присутствуют в обеих МНСЧБ, но количество нулей в первой 
матрице несоизмеримо превосходит количество нулевых значений второй матрицы, 
а элементы, равные 5, 6, 7, 8, в первой матрице отсутствуют вообще. 
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Обнаружение фальсификации 


Будем считать, что в нашем распоряжении имеются фотоснимки, полученные 
современными цифровыми фотокамерами, то есть каждый из них — это полностью 
восстановленное после ЛРЕС-сжатия изображение. Тогда соответствующие МНСЧБ 
будут иметь вид, аналогичный представленному на рис. 2(6). 

Пусть часть фотоснимка, который будем называть основным изображением 
(ОИ), заменяется частью ЦИ, далее называемой вклейкой, или замещающей 
областью (ЗО), хранящегося в формате без потерь, причем для большей наглядности 
получаемых ниже выводов никакая последующая обработка изображения не 
производится. Такое фальсифицированное изображение (фотомонтаж) сохраняется 
без потерь. Пример на основе изображения САМЕКАМАМ.7Т\ТЕ, намеренно явно 
демонстрирующий вклейку, представлен на рис. 3З(а). При построении МНСЧБ 
полученного фотомонтажа части, отвечающие ОИ и ЗО, будут отличаться по 
количеству нулевых СЧ блоков: подобласть МНСЧБ, которая соответствует вклейке, 
будет содержать большое количество нулей и значений, мало отличающихся от нуля 
(на рис. 3(6) эта часть ограничена). 

Таким образом, для детектирования ЗО необходимо выделить в МНСЧБ связные 
области, содержащие нули и близкие к нулю значения. Однако такие совокупности могут 
отвечать: а) вклейкам; 6) областям ОИ, где малое количество нулевых СЧ обязано 
округлениям, проводимым после частичного восстановления изображения (их наименьшие 
СЧ будут незначительно отличаться от нуля, что не свойственно вклейке); в) областям ОИ, 
у которых после квантования и ЧВ уже имелись малое количество нулевых СЧ, то есть 
областям, содержащим значительное количество контуров. На рис. 3(6) примеры 
областей 6) и в) ограничены другим цветом. Очевидно, что различить области а) и 6), 
сравнивая сингулярные спектры блоков, входящих в них, не представляет труда (табл. 4 — 
«вклейка» и «основное изображение» области 1, 2). Однако поведение СЧ блоков 3-ей 
части ОИ качественно неотличимо от поведения СЧ блоков ЗО и никак не выделяет 


вклейку. 


Рисунок 3. — Фальсифицированное изображение (а) и соответствующая МНСЧБ (6) 
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Таблица 4 — Сингулярные спектры различных областей фальсифицированного 


изображения 
Блок Сингулярные числа 
изображения 
(21,19) 939,0048 106,2395 82,7256 63,2449 26,6378 17,8791 4,1803 0,1620 
8 (22,20) 819,4101 214,7784 199,6390 85,0353 45,4888 29,4758 21,3788 6,2910 
я (23,20) 446,5187 227,7229 72,7133 53,4687 16,6807 4,3490 3,1316 1,1126 
5 (23,22) 1095,5 197,5 120,5 94,9 28,1 19,8 19,2 1,3 
г (24,23) 659,2959 86,3697 16,8607 10,6777 6,6923 5,8589 2,7147 0,6267 
(25,22) 1099,9 289,7 134,4 72,8 50,5 33,3 8,8 2,2 
(20,3) 1175,8 96,4 57,9 21,8 7,2 0,7 0,4 0,4 
И (20,4) 1122,3 106,1 63,1 38,6 14,9 0,7 0,6 0,1 
Е 1 (21,2) 702,7781 28,2405 19,2169 13,1577 0,7773 0,5476 0,2373 0,0735 
я (21,3) 675,3140 46,6384 39,9138 16,1638 8,7743 0,4852 0,4037 0,2509 
8. (22,3) 1035,5 111,0 3,.2 22,4 1,8 0,8 0,3 0,1 
8 (20,25) 1145,0 98,4 79.2 32,7 19,0 0,5 0,3 0,1 
ь 2 (21,26) 1009,1 181,5 67,2 24,9 15,8 0,9 0,5 0,3 
8 (22,27) 768,9128 195,9938 95,1066 79,8180 13,2214 0,6426 _ 0,4150 0,1462 
а (11,17) 1001,6 2237 180,4 1367 106,6 60,3 6,6 0,5 
© (11,18) 651,0872 124,1137 90,5530 59,1854 30,2031 21,4353 2,0865 0,6501 
3 (12,17) 1124,2 244,3 144,9 137,5 91,0 52,4 4,0 1,3 
(12,19) 1026,1 220,8 151,9 117,3 59,8 13,7 7,2 1,4 


Проведем /РЕС-сжатие для смонтированного изображения, но элементы мас- 
сива нормализации возьмем меньше, чем использованные при первом квантовании 
коэффициентов ОСТ ОИ. Из сказанного выше вытекает, что второе квантование и 
ЧВ не должны качественно изменить картину для сингулярных спектров блоков 3-ей 
области ОИ. Действительно, отсутствие нулевых СЧ в блоках этой части, 
значительное отличие последних по порядку СЧ от нуля говорит о том, что эти 
блоки содержат большое количество контуров, то есть значительными являются 
коэффициенты ОСТ, соответствующие высоким и средним частотам. Если эти 
коэффициенты не обнулились при первом квантовании, то болышая их часть 
останется ненулевыми и при втором, где коэффициенты квантования будут меньше. 
Для вклейки процесс квантования коэффициентов ОСТ будет первым, поэтому 
проявится картина результатов первого квантования — значительное увеличение 
количества нулевых СЧ в ее блоках. Проиллюстрируем вышесказанное на рас- 
сматриваемом примере. 


Рисунок 4 — МНСЧБ частично восстановленного после /РЕС-сжатия 
смонтированного изображения 
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Рисунок 5 — Фальсифицированное изображение (а), область МНСЧБ, отвечающая ЗО 
до сжатия (6) и область МНСЧБ, отвечающая ЗО после ЧВ (в) 


При ЧВ смонтированного изображения, как и ожидалось, для области 3 ОИ кар- 
тина количества нулевых СЧ в блоках изменилась незначительно (нулевые элементы 
матрицы стали | или 2). Для вклейки нулевые значения МНСЧБ изменились в 
диапазоне от | до 5 (рис. 4). Заметим, что в предложенном примере вклейка сама по 
себе содержит много контуров. Если же вклейка представляет из себя часть гладкого 
изображения с малым количеством или отсутствием контуров (рис. 5(а)), то 
отличить ее от области типа 3 ОИ будет значительно легче: элементы МНСЧБ будут 
больше отличаться от нуля после ЧВ изображения (рис. 5(в)). 

Таким образом, если ОИ является результатом полного восстановления 
изображения после ГРЕС-сжатия, а ЗО — часть изображения, хранившегося без потерь, 
причем фальсифицированное изображение сохранено без потерь (ТЕ, ВМР и т.д.), то 
основные шаги предлагаемого метода определения наличия вклейки и локализации ее 
области будут следующими: 

1. Матрица анализируемого ЦИ разбивается стандартным образом на блоки 
8х8, для каждого блока находится сингулярный спектр, результатом чего является 
МНСЧБ М.. 

2. В М, выделяются связные области О,О.,...О»„ (не обязательно прямо- 


угольные), абсолютное большинство элементов которых имеют нулевое значение. Эти 
области могут соответствовать вклейкам. 
3. Из О,,О.,....О„ удаляются те, для которых значительная часть минимальных 


СЧ блоков мало отличается от нуля. Пусть О, ,О,,,..,О, — оставшиеся области. Если 
Р 


множество оставшихся областей пусто, то вклейки отсутствуют. 

4. Анализируемое ЦИ подвергается ЛРЕО-сжатию, при котором элементы 
массива нормализации будут меньше, чем для /РЕС-сжатия ОИ, с последующим ЧВ. 
Для матрицы ЧВ-изображения получаем МНСЧБ М.. 


5. В матрице М, выделяем подобласти Охк,Оь,..,Ок,, отвечающие 
Ио О, в М,. Вклейкам соответствуют те из Ок,Оь......Ох, ‚ которые содержат 


значительное число элементов из множества {3, 4, 5, 6, 7, 8}. Если такие области не 
обнаруживаются, то вклейки отсутствуют. 

Замечание. При формировании монтажа ОИ не использовались дополнительно 
операции обработки изображения для маскировки вклейки. Однако их применение 
не изменит характерных особенностей СЧ, выявленных выше, отличающих блоки 
ОИ от блоков ЗО, ведь основа этих отличий лежит в том, что ОИ было сохранено с 
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потерями, а вклейка без потерь. Заметим лишь, что размывка контура ЗО может 
привести к изменению картины полностью восстановленного после ЛРЕО-сжатия 
изображения в областях ОИ, прилегающих к вклейке, что может помешать точному 
определению ее границы. 


Выводы 


Настоящая работа является первым шагом на пути формирования нового подхода 
к проблеме детектирования фальсификации цифрового изображения, использующего в 
качестве математических инструментов теорию возмущений и матричный анализ, что 
никогда не делалось ранее. Предложенный в работе метод на основе анализа 
возмущений сингулярных чисел блоков матрицы ЦИ дает возможность не только 
установить сам факт несанкционированного изменения цифровой фотографии, но и 
локализовать область такого изменения, не привлекая для этого никакие дополнительные 
технические средства, а также средства обработки полученного изображения. Вся анали- 
зируемая информация предоставляется непосредственно матрицей ЦИ. 

Хотя предложенный метод рассматривает пока лишь частный случай фальси- 
фикации ЦИ, однако универсальность используемого подхода позволяет надеяться на 
возможность уйти в дальнейшем от конкретики фальсификации, анализируя лишь 
изменения СЧ матрицы изображения в результате возмущающего воздействия. 
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А.А.Кобозева 

Використання особливостей збурень сингулярних чисел матрищ числового зображення для 
виявлення його фальсифкаци 

Пропонуеться новий шдхд до розв’язання проблеми виявлення фальсиф1кацй цифрового зображення, 
який базуеться на анал!з1 збурень сингулярних чисел блоюв матриц! зображення. Новий шдх!д дае 
можлив1сть локалзовати область несанкцоновано! змни поданного зображення. 


А.А. КоБозеуа 

АррНсайоп о? Ппазе Майлх Эшошаг Уашез О15игапсе$ гог Ппасе Еогоегу Ваесйоп 

Меу\ арргоасВ ® а ргоет оЁ пптазе Ююгсегу деесйоп Базе оп апа]у$1$ оЁ зшечаг уаез 41иФапсез оЁ 
плазе таблх 15 ргорозед ш 1$ агае. ТЬ1$ арргоасВ епаЫез {0 |осайе гапое оЁР ипалопяей уапаНоп оЁ 
опеша! ппазе. 
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